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Abstract 

The possibility of exciting the g-modes of a solar-type star as a consequence 

of the gravitational interaction with a close companion (a planet or a brown 

dwarf) is studied by a perturbative approach. The amplitude of the emitted 

gravitational wave is computed and compared with the quadrupole emission 

of the system, showing that in some cases it can be considerably larger. The 

effects of radiation reaction are considered to evaluate the timescale of the 

emission process, and a Roche lobe analysis is used to establish the region 

where the companion can orbit without being disrupted by tidal interactions 

with the star. 
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I. INTRODUCTION 



In recent years a significant number of systems composed by a solar-type star and one 
or more orbiting companions has been discovered in our neighbourhood. According to the 
Extrasolar Planets Catalog ffl], at the time we write this paper 58 extrasolar planetary 
systems have been confirmed, and this number is destinated to grow in the near future (see 
0] for a recent review). 46 of the observed orbiting objects are planets with masses ranging 
between [0.16-ll]Mj, where Mj is Jupiter's mass; 12 are bigger, up to 60 Mj, and include 
super-planets and brown dwarfs. In the following, we shall indicate all the orbiting objects 
as 'planets'. 

The newly discovered systems exhibit some unexpected features; 22 planets move on 
elliptic orbits with eccentricities e > 0.3, larger than the largest eccentricity encountered in 
the solar system (Mercury, e = 0.2), and many of them are incredibly close to the central 
star, much closer than expected. For instance the planet orbiting 51 Pegasi, discovered 
in 1995 0, ||, has a period P = 4.23 days and orbits its star at a distance a = 0.05 
AU. More than two thirds of the planets discovered up to now are orbiting their host star 
much closer than Mercury orbits the Sun (P = 88 days, a = 0.39 AU), and even closer 
than 51 Pegasi's planet. These orbital properties challenge the existing theories on planet 
formation and evolution, and suggest an interesting possibility: there may exist planets on 
orbital radii sufficiently small to excite some g-modes of the star by tidal interaction. This 
issue has been discussed by Terquem et al. f| in connection with the short orbital period 
of 51 Pegasi's planet. They calculated the dynamical tides raised by the planet on the 
star, and the dissipation timescales due to turbulent viscosity in the convective zone and 
to radiative damping in the radiative core. In this paper we shall approach the problem 
of the excitation of the g-modes from a different point of view: by using a perturbative 
approach, justified by the fact that the mass of the central star is much bigger than that 
of the planet, we shall compute the gravitational signal emitted by a system in which the 
planet moves on a circular orbit of radius R , close to that which would correspond to the 
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resonant excitation of a g-mode, R res - We shall extend the work presented in a previous 
paper ( 0, to be referred to hereafter as Paper I), were we have studied the quadrupole 
emission due to the orbital motion of the observed planetary systems, and started to explore 
the possibility for a planet to excite some stellar modes. First of all, in section II we shall 
verify by a Roche-lobe analysis whether the planet can move on the orbit R = Rq without 
being disrupted by tidal interactions. In section III, using the quadrupole formalism we 
shall compute the characteristic amplitude of the gravitational signal, Iiq, emitted by the 
system because of the orbital motion of the planet on Rq. As usual, the two bodies will be 
treated as pointlike masses, and the effects of the planet on the internal structure of the 
star will be neglected. These effects will be considered in section IV, where we shall assume 
that the planet induces a small perturbation on the gravitational field of the star and on its 
thermodynamical structure. We shall integrate the perturbed Einstein's equations coupled 
to the hydrodynamical equations, having the stress energy tensor of the planet as a source, 
and compute the characteristic amplitude of the emitted gravitational signal, /i^, which 
will be compared with hq. In this way we will be able to evaluate the enhancement in the 
radiation emitted by the system as a consequence of the g-modes excitation. In Section V 
we shall discuss how long can a planet move on an orbit close to resonance before radiation 
reaction effects lead the planet off resonance. Conclusions will be drawn in section VI. 

II. WHICH STELLAR MODES CAN BE EXCITED BY AN ORBITING PLANET 

In this section we shall verify, by a Roche-lobe analysis, if a planet can orbit on a radius 



such that the keplerian orbital frequency, uj^ = yG(M+ + M p )/Rq, is half the frequency of 

a given mode of the star uji. As discussed in [|J, [[?] and ||, in this case the mode would 
be excited by the dynamical tides raised by the planet. We shall consider as a model a 
polytropic star, p = Kt l+1 l n with n = 3, e /po — 5.53 • 10 5 and adiabatic exponent 
7 = 5/3. Choosing the central density eo = 76 g/cm 3 , this model gives a star with the 
same mass and radius as the Sun. In Paper I we showed that if one considers a polytropic 
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star with n = 2 and mass equal to that of the Sun, the p-modes and the fundamental mode 
cannot be excited by a planet, because it would be tidally disrupted. The same holds for the 
more appropriate, though still simplified, n = 3 model considered in this paper, therefore 
we shall deal only with the excitation of the g-modes. 

We shall consider three companions for the central star: two planets with the mass of the 
Earth and of Jupiter, and a brown dwarf of 40 jovian masses, Me, Mj and Mbd, respectively 
Smaller planets produce gravitational signals that are too small to be interesting. Let us 
assume that a planet orbits the star at a constant orbital radius R l Tes , in a position to excite 
the mode gi. Following the procedure described in Paper I it is possible to determine, by a 
Roche lobe analysis, what is the maximum radius the planet can have in order not to overflow 
its Roche lobe. Since we have fixed the mass of the planet, this constraint also determines 
the minimum value of the planet's mean density, p m i n , compatible with the excitation of 
that mode. The results of this analysis are shown in table 1, where we tabulate the value of 
Pmin, expressed in units of the mean density of the Sun, for the three considered companions. 
It should be stressed that in the Roche lobe analysis the ratio p m in/ P*, for each given mode, 
depends only on the ratio between the mass of the planet and that of the star. From table 
| we see that, since PEarthl P& = 3.9 and pj/p® = 0.9, a planet like the Earth can orbit 
sufficiently close to the star to excite g-modes of order higher or equal to n = 4, whereas a 
planet like Jupiter can excite only the mode gio or higher. 

According to the brown dwarf model ( "model G" ) by Burrows and Liebert an evolved, 
40 Mj brown dwarf has a radius Red — 5.9 • 10 4 km, and a corresponding mean density 
Pbd = 88 g-cm~ 3 ; consequently Pbd/Pq = 64, and this value is high enough to allow a 
brown dwarf companion to excite all the g-modes of the central star. However, we also need 
to take into account the destabilizing mechanism of mass accretion from the central star. 
This imposes a further constraint, and this is the reason why the slots corresponding to the 
excitation of the g-modes lower that g4 in the last row of table | are empty. 
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III. THE QUADRUPOLE EMISSION OF A BINARY SYSTEM 



In this section we shall use the quadrupole formalism to compute the gravitational signal 
emitted by a binary system. We shall assume that a planet orbits its host star in circular 
motion, and that the radius of the orbit corresponds to the excitation of one of the g-modes 
allowed by the Roche lobe analysis discussed in section II. As usual, both the star and the 
planet will be treated as pointlike masses, with no reference to their internal structure. From 
the reduced quadrupole moment 

Q kl = // (x k X l - |X| 2 V (3.1) 
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where /i = MpM*/ (M p + M*) and X = x* — x p is the relative position vector, the non 
vanishing TT-components of the emitted wave can easily be computed 

hj^(t, x) = H [p ik p 3l - Q kl (r)\ T=t ^_ (3.2) 

where Pik = 5 l k —riink is the projector onto the 2-sphere r = const, and n is the radial unit 
vector. It should be stressed that eq. ( |3.2| ) allows to compute the radiative contribution due 
solely to the orbital motion of the system. For a circular orbit the wave is emitted at twice the 
keplerian frequency |[L0|| , and from the two independent wave components h±(2u>k,r, , }),ip) 



we compute the characteristic amplitude ]TTJ 



12 1 



|2 , ir |2 



h Q (2u k ,r) = \l-J— /ain \h+\* + \h-\* =^\j-{u k Ro)\ (3.3) 



M v A . 2 , D , a 



r V 15 



where the factor y 2/3 takes into account the average over orientation, and f2 is the solid 
angle. This amplitude will be compared with that determined by the perturbative approach. 
In table |Tj we give the values of Hq for the three companions orbiting the sun-like star which 
we use as a model. The orbital radius is chosen to correspond to the excitation of a permitted 
g-mode (cfr. table |). The planetary systems are assumed to be located at a distance of 10 
pc from Earth. In the last row we give the frequency (in /iHz) of the g-modes, which is also 
the frequency at which waves are emitted. 
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IV. THE GRAVITATIONAL SIGNAL EMITTED BY A PERTURBED STAR 



In this section we shall briefly introduce the equations of general relativity which describe 
the perturbations of a star excited by a planet moving on a circular orbit. Unlike the 
quadrupole formalism discussed in section III, the perturbative approach takes into account 
the internal structure of the star and the modification induced on it by the interaction with 
the planet, and it allows to evidentiate the contribution to the emitted wave due to the 
excitation of the stellar modes. When the orbit is circular, only the polar, nonaxisymmetric 
perturbations are excited, and the metric appropriate for their description can be written 



in the following form [I2| 
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where the functions [Ni m (cj, r), T; m (cj, r), Vi m (u, r), Li m (cj, r)) describe the radial part of the 
perturbed metric, and VJ m (i?, <p) are the scalar spherical harmonics. The metric functions 
v{r) and ^{f) describe the unperturbed spacetime and are determined by numerically inte- 
grating the TOV equations for hydrostatic equilibrium for the chosen equation of state. The 
functions [N, T, V, L] have to be found by solving the perturbed Einstein equations outside 
the star and inside, where they couple to the hydrodynamical equations. Inside the star, 
after separating the variables and assuming a nonbarotropic equation of state, the perturbed 
equations can be reduced to the following set, from which the hydrodynamical variables have 
been eliminated 
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In eqs. (|4.2j ) the harmonic indices have been suppressed, we have replaced the 

function V by X = nV, with n = (I — + 2)/2, and the function T by 



G = u, 
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We have also defined Q = (e + p) / 7p, where 7 is the adiabatic exponent 

(e + p) ( dp\ 
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P 



de J 



entropy=const 

and e and p are the energy density and the pressure of the fluid composing the star. Since 
in this paper we apply the relativistic theory to nonbarotropic stars, we choose, as in ref. 
|i~3|l , 7 = const = 5/3. The system of eqs. ( [4.2|) can be solved numerically by integrating 
the two independent solutions which satisfy the regularity condition at the center, and 
superimposing them in such a way that the perturbation of the pressure vanishes at the 



boundary, as discussed in [14 



Outside the star, e and p vanish, and a source term given by the stress-energy tensor 
of the moving planet must be added on the right-hand-side of the Einstein equations. Eqs. 



( OP reduce to the Zerilli equation (Tj| for the Schwarzschild perturbations 

\ d 2 2 2(r -2Mjn 2 (n + l)r 3 + 3M 2 r 2 +9M 2 nr + 9M*]) „ . . n , , 
{^ + » r , {nr + 3MJ2 ^Z lm {^r) = S lm {u,r), 



(4.4) 



where r* = J r e U+I12 dr, and M* is the mass of the star. The Zerilli function matches 
continuously with the solution of eqs. (|4.2|) in the interior, i.e., at the surface of the star 



Zim(w,R) = p (—X lm (u,R)-RL lm (u,R)), (4.5) 

and similarly for its first derivative (see ref. [JTJ] for details). The source term, Si m (uj,r), is 
derived from the stress-energy tensor of the planet, considered as a pointlike particle moving 
on a circular orbit around the star in the equatorial plane 

From the geodesic equations ^ = E (l — , where the planet's energy per unit rest 



E={1 .m( l -m- 1 ' 2 . ( 4.7) 



mass is 
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In terms of these quantities the source term can be written as 

[ 87r(r - 2i\4)(12M 2 + r 2 n 2 + 3nrM, - 6rMj (0) 
Si m {w,r) = M p < — (4.8) 
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In order to find the amplitude of the emitted wave we follow the same approach used by 
Kojima ||16|| , who first calculated the gravitational signal emitted by a particle in circular 
motion around a star. In particular, he considered the excitation of the fundamental mode 
of a neutron star by an orbiting particle, showing that a sharp resonance occurs if the 
frequency of the f-mode is twice the orbital frequency, and that the characteristic wave 
amplitude emitted at that frequency can be up to 100 times larger than that evaluated by 
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the quadrupole formula. As discussed in section II, if the central star is of solar-type only 
the g- modes can be excited by a planet, and therefere we shall focus on these modes. 

The solution of eq. (|4.4|) can be constructed by the Green-function technique as follows. 
We first integrate the interior equations fl4.2|) from r = to r = R, where we compute 
the Zerilli function Zi m (R) as given in eq. (f4.5|), and its first derivative. We then 



construct, by numerical integration, a solution of the homogenous Zerilli equation, Z} m1 
which satisfies the boundary condition Zl m (R) = Zi m (R); at radial infinity Z 1 behaves 
as a superposition of ingoing and outgoing waves, Z\ m ~ A\ m e~ lujr * + Af^e +lujr * . Let us 
consider a second solution of the homogenous Zerilli equation, Zf m , which behaves as a pure 
outgoing wave at infinity, i.e. Zf m ~ e +lujr *. The general solution of the non- homogeneous 
Zerilli equation ( f4.4|) , which satisfies the physical requirement of pure outgoing radiation at 
infinity, Zi m {r — > oo) ~ Zf^{ijj)e +lwr * , and the matching condition at the surface, can be 
written in terms of these two functions 



1 r r r * 

J lm 



w lm 

where Wi m is the wronskian of the two solutions 



{ Z lmJ R Z lm S lm dr *- Z lm J Z? m S lm dr^ (4.10) 



Wim = {zl m Zf mr ^ — Zf m Zl mr ^ — 2iujA\ m . (4-11) 



The amplitude of the Zerilli function at infinity therefore is: 



I roo 

^f(") = / ZLShndr*. (4.12) 



2iuA\ m J R 

Because of the form of the source term ( |4.8|) , (|4.9|) , this amplitude can be written as 

ZZ l = M p ZZ\mw kl R )S(u - mcu k ), 

It should be noted that if we restrict our analysis to I = 2, the source is non zero only 
for m = ±2. Finally, the amplitude of the outgoing gravitational radiation at infinity can 
be computed by using the relation between Z out and the two wave components h± in the 
radiative gauge |15| as follows 
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h + (2u k ,r) + ih_(2uj k ,r) 



(4.13) 
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where is the s = 2 spin-weighted spherical harmonic. As in section III, we shall define 
the relativistic characteristic amplitude, h 1 ^ 2 , to be compared with fiQ, by 

_ Mp 2^n(n + l)\2^(2uj k ,R ) 
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where we have used the ort honor mality of the spin-weighted harmonics, and, in summing 



ryout 
^22 



yout 
^2,-2 



over to, the symmetry property 

We have computed /i^jr 2 assuming that the planet moves on a circular orbit of radius 
Rq = R l res + AR, where R l res is the orbit corresponding to the excitation of the mode for 
the modes allowed by the Roche lobe analysis. We find that, as the planet approaches the 
resonant orbit, h 1 ^ 2 grows very sharply. It is instructive to plot the ratio h l ^ 2 (R )/hq(R l res ) 
as a function of AR, to see how much the amplitude of the emitted wave grows with respect 
to the quadrupole emission, because of the excitation of a g-mode. In figure 1 the logarithm 
of this ratio is plotted for the modes g 4 , g 7 , and g 10 as a function of logAi?, showing a 
power-law behaviour nearly independent of the order of the mode. It should be stressed that 
h l ft 2 (R ) I l hQ(R res ) is independent of the mass of the planet, but depends, of course, on the 
selected stellar model. From figure 1 we deduce that, in principle, as the planet approaches 
a resonant orbit the amplitude of the emitted wave may become significantly higher than 
that emitted because of the orbital motion {Jiq). Thus a relevant question to answer is how 
long can a planet move on an orbit close to a resonance, before radiation reaction effects 
move it off resonance. This issue will be discussed in the following section. 



V. THE EFFECT OF RADIATION REACTION 

The loss of energy in gravitational waves causes a shrinking of the orbit of a planetary 
system, and the efficiency of this process increases as the planet approaches a resonant 
orbit. We shall now compute the time a planet takes to move from an orbit of radius Ro = 

10 



R\ es + AR, where the amplification factor KrIKq has some assigned value, to the resonant 
orbit R % res i because of radiation reaction effects. This timescale will indicate whether a planet 
can stay in the resonant region long enough to be possibly observed. On the assumption 
that the timescale over which the orbital radius evolves is much longer than the orbital 
period (adiabatic approximation), the orbital shrinking can be computed from the energy 
conservation law 

dEcw 



dE\ 
( ~dt } + 



dt 



0. 



(5.1) 



where E is the energy per unit mass of the planet as given in (|4.7|) , and ^ dE j t w ) is t ne energy 



emitted in gravitational waves, which can be computed in terms of the wave amplitude Q4.12p 
as follows (cfr. e.g. JT7[) 
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Since (f ) = (^) / (^), using eq. Q and eq. Q, eq. O gives 



dRn 



2Rl (1-3MJR ) 3 / 2 ldE GW 



\ dt / M p M+ (1 - QMJR ) \ dt / ' 
from which the time needed for the planet to reach the resonant orbit can be computed 

(1 - 6MjRa) dR 



(5.3) 



AT 



M P M+ r w r t 



(5.4) 



2 Jr^+ar {^-) (1 - 3MJR ) 3 / 2 Rl 

It should be noted that since (^g*) is proportional to M^, AT is longer for smaller 
planets. We have computed AT for the three companions Mg, Mj and Mbd, by the following 
steps: 

- For each g-mode allowed by the Roche lobe analysis we find the radius of the resonant 
orbit Rl es . 

- We assume that each companion orbits the star at a distance Ro = R l res + AR such that 
the amplification factor A = hu(R )/hQ(R l res ) has an assigned value, and compute the 
corresponding energy radiated in gravitational waves, (^ dE ^ t w 
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- We compute AT, which tells us how long can the companion orbit the star in the resonant 
region between R and R l res , emitting a wave of amplitude higher than A hq(R l res ). 
The results are summarized in table fTJ. These data have to be used together with those in 
table y as follows. Consider for instance a planet like the Earth, orbiting its star on an orbit 
resonant with the mode g 4 . According to the quadrupole formalism, which does not consider 
changes in the internal structure of the star and therefore does not include the resonant 
contributions to the emitted radiation, it would emit a signal of amplitude fiQ(g±) = 3-1CT 26 , 
at a frequency vow = 1.1 • 10~ 4 Hz (table U first and last row, respectively). The data of 
table |n|, which include the resonant contribution, indicate that before reaching the resonant 
orbit Rf es , the Earth-like planet would orbit in a region of thickness AR = 8.7 km slowly 
spiralling in, emitting waves with amplitude Hr > lO/ig^) = 3- 10 -25 , for a time interval of 
8.7- 10 6 years, and that while spanning the smaller radial region AR = 1.3 km, the emitted 
wave would reach an amplitude Hr > 50/iq((?4) = 1.5 • 10~ 24 , for a time interval of 4.7 • 10 4 
years. 

A jovian planet, on the other hand, could only excite modes of order n = 10 or higher, 
which would correspond to a resonant frequency of vqw — 5-7 • 10~ 5 Hz and a gravitational 
wave amplitude greater than 6 • 10 -23 for 520 years (AR = 42 m), or greater than 3 • 10~ 22 
for ~ 3 years (AR = 6m). From these data we see that the higher the order of the mode, 
the more difficult it is to excite it, because the region where the resonant effects become 
significant gets narrower and the planet transits through it for a shorter time. 

Much more interesting are the data for a brown dwarf companion. In this case, for 
instance, the region which would correspond to the resonant excitation of the mode g4 
(AR = 8.7 km), with a wave amplitude greater than 3.9 • 10~ 21 , would be spanned in 690 
years, whereas the emitted wave would have an amplitude greater than ~ 2- 10 -20 (AR = 1.3 
km) over a time interval of ~ 4 years. 
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VI. CONCLUDING REMARKS 



The study carried out in this paper has been motivated by the recent discovery of sev- 
eral extrasolar planetary systems very close to Earth, with companions orbiting the host 
star much closer than previously expected. This discovery suggests the possibility that a 
companion may be found orbiting at such close distance to excite one of the g-modes of 
the star, and a Roche-lobe analysis indicates that this may be possible, especially if the 
companion is a brown dwarf. The frequency of a g-mode in a sun-like star is of the order 
of a few digits in 10~ 4 Hz, and would be in the bandwith of the spaceborn interferometer 
LISA. Thus, our aim was to understand how much such resonant process may enhance the 
emission of gravitational waves of the system with respect to the radiation it emits because 
of the time varying quadrupole moment associated to the orbital motion. This enhance- 
ment has been evaluated by integrating the equations that describe the perturbation the 
companion induces on the star, and computing the amplitude of the emitted wave when the 
planet moves close to an orbit corresponding to the resonant excitation of a g-mode. The 
amplitude increases sharply as the orbit approaches the resonant one; however, in order to 
get a chance to observe one of these signals emitted by a system in our vicinity, it is also 
needed that it stays above a certain threshold long enough. For this reason we have also 
computed the effects of radiation reaction on the shrinking of the orbit and the associated 
timescales. The main information we extract from our study is that only close brown dwarf 
companions may produce signals of some relevance; a brown dwarf of 40 Jovian masses, 
orbiting a star located at a distance of 10 pc from Earth, if close to an orbit resonant with 
the mode g 4 , may emit radiation of amplitude greater than ~ 2 • 10~ 20 at a frequency of the 
order of ~ 10~ 4 Hz, spanning a region AR =1.3 km over a time interval of ~ 4 years. Of 
course the amplitude of the wave increases with the mass of the brown dwarf and decreases 
with the distance, and it will be interesting to see whether astronomical observations will 
identify any such system in our vicinity in the future. 

As a last remark, we would like to mention that things would not be significantly different 
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had we considered more realistic models of solar-type stars, such as those discussed in fTE] 



fj~9|| ; they predict oscillation frequencies that are a little higher than those of our simple 
polytropic model; however, our frequencies are at least qualitatively correct to give order- 
of-magnitude estimates for the effects of resonant excitation of the stellar modes and for the 
excitation timescales. 
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TABLES 

TABLE I. The minimum mean density, p m in/PQ, that a planet should have in order to be 
allowed to move on an orbit which corresponds to the excitation of a g-mode, without being 
disrupted by tidal interactions. The values of Pmin/p& are tabulated for three planets with mass 
equal to that of the Earth (Me), of Jupiter (Mj) and of a brown dwarf with Mbd = 40 Mj. 
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3.68 
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2.08 


1.65 


1.34 
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0.93 



TABLE II. The amplitude of the gravitational signal emitted when the three companions 
considered in table 1 move on a circular orbit of radius R l res , such that the condition of resonant 
excitation of a g-mode is satisfied, is computed by the quadrupole formalism for the modes allowed 
by the Roche-lobe analysis. The planetary systems are assumed to be at a distance of 10 pc from 
Earth. In the last line the emission frequencies vgw are given. 
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TABLE III. In column 2 we give the frequency, vgw-, of the wave emitted when a companion 
moves around the host star on an orbit resonant with a g-mode allowed by the Roche lobe analysis 
(see text). Because of resonant effects, the amplitude of the wave is amplified by a factor greater 
than A (column 3) when the companion spans a radial region of thickness AR (column 4), which 
is the same for all planets. In the last three columns we give the time interval AT needed for the 
three companions to span the region AR and reach the resonance, because of radiation reaction 
effects. 
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FIGURES 




AR (km) 

FIG. 1. The logarithm of the ratio h l ^ 2 (R % res + AR)/hQ(R l res ) is plotted as a function of AR 
for the modes g4, g7, and gio (see text). 



19 



